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ABSTRACT 

We show that the significantly different effective temperatures (T e g) achieved by the luminous 
blue variable AG Carinae during the consecutive visual minima of 1985-1990 (T c ff — 22,800 K) and 
2000-2001 (T c s ~ 17,000 K) place the star on different sides of the bistability limit, which occurs in 
line-driven stellar winds around T e g ~ 21, 000 K. Decisive evidence is provided by huge changes in the 
optical depth of the Lyman continuum in the inner wind as T e g changes during the S Dor cycle. These 
changes cause different Fe ionization structures in the inner wind. The bistability mechanism is also 
related to the different wind parameters during visual minima: the wind terminal velocity was 2-3 
times higher and the mass-loss rate roughly two times smaller in 1985-1990 than in 2000-2003. We 
obtain a projected rotational velocity of 220 ± 50km s _1 during 1985-1990 which, combined with the 
high luminosity (L* = 1.5 X 10 6 -Lq), puts AG Car extremely close to the Eddington limit modified by 
rotation (ilT limit): for an inclination angle of 90°, Fq ^ 1-0 for M < 60 Mq. Based on evolutionary 
models and mass budget, we obtain an initial mass of ~ 100 Mq and a current mass of ~ 60 — 70 Mq 
for AG Car. Therefore, AG Car is close to, if not at, the fir limit during visual minimum. Assuming 
M = 70 Mq, we find that Tq decreases from 0.93 to 0.72 as AG Car expands toward visual maximum, 
suggesting that the star is not above the Eddington limit during maximum phases. 
Subject headings: stars: atmospheres — stars: mass loss — stars: variables: other — supergiants — 
stars: individual (AG Carinae) — stars: rotation 



1. INTRODUCTION 

Albeit rare, massive stars are the main contributors 
to the input of ionizing photons, energy, and momen- 
tum into the interstellar medium, and are responsible 
for a significant fraction of the chemical enrichment of 
their host galaxy. Massive stars evolve on relatively 
short timescales (a few 10 6 years), and their evolution 
is strongl y influenced by their stron g mass loss and rota- 
tion (e.g. IMevnet fe Maedeill2005al) . 

Tremendous advancement in the understanding of the 
evolution of massive stars has been achieved in recent 
decades (lContilll976t iSchaller etallll992j: IMevnet et all 
119941: iMaeder fe Mevnetl 119941 iLanger et al.l I1994D . ~hi 
particular, we now have insights into the effects of rota- 
tion on the evolution of massive stars (jMaeder fe Mevnetl 
I2000bl) . The models predict the existence of a short- 
lived, transitional stage, usually referre d to as the 
Luminous Blue Variab l e (LB V) phase (jContil 119841: 
iHumphrevs fe Davidson! I1994D . during which the star 
has a high mass-loss rate (M ~ 10~ 5 — 10 -3 M Q yr _1 ). 
In the current picture of stellar evolution, LBVs are 
rapidly evolving massive stars in the transitory phase 
from being an O-type star burning hydrogen in its core 
to a hydrogen-poor Woff-Rayet, helium-core burnin; 
star (IHumphrevs fe Davidson! | 1994t IMaeder fe Mevne 
I2000bl: IMevnet fe Maederll2000[ 120031) . The LBV phase 
is likely an unavoidable evolutionary stage which mas- 
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sive stars with initial mass Af; n i > 40 Mq will ex- 
perience during their lifetime before becoming WRs 
and e xploding as a supernova ([Humphr eys fe Davidson] 
1994). Whether LBV progenitors have an upper limit 
for Mj n j, as suggested by th e latest evolutionary mod- 
els ([Mevnet fe Maedeil 120031 ). still remains to be seen 
Several LBVs, such as E ta Car ([Davidso n fe Humphreys! 
[19971 IHillier et al.1 12001 . AG C ar (IGroh et al. M2009bD . 
and th e Pistol star (Figer etaTI I1998D. among oth- 
ers (seelHumphrevslll978t IHumphrevs fe Davidsonlll994l ; 
Ivan Gendered 120011; IClark et al.l I2005D . are known to 
have bolometric luminosities in excess of 10 6 Lq, which 
suggests a progenitor with Mi n i > 100 Mq. 

Following recent determinations that the mass- 
loss rates of O-type stars are a factor of 3-20 
lower than t hose adopted by the e volutionary models 
cited above dCrowther et al.l 120021; I Hillier et al.l 
Massa et al.ll2003l; Kramer et al.ll200a|Evans et al.l 



2003 



2004; 



Bouret et al.l 120031 120051; iPuls et al.l I2006D . and the 

detection of massive nebulae around LBVs, there is 
growing evidence supporting the idea that the evo- 
lution of massive stars is even more dominated by 
strong m ass loss during the LBV phase than previously 
thought ([Smith fe Owockil 12006). Surprisingly, some re- 
cent work sugge sts that some core-collapse SNe have 
LBV progenitors (IKqtak fe Vinkll2T)0llSmith et al.ll2007t 
IGal-Yam et "all I2007t iTrundle et al.l (20081 ). which dra- 
matically enhances the cosmological importance of LBVs 
and poses a great challenge to the current paradigm 
of massive star evolution. Some of these type-II SNe 
occurred just a few years after the progenitor suffered 
a giant mass ejec tion a la Eta Carinae in the 1840's 
(Smit h et al.ll2007[) . raising the following provocative hy- 
pothesis as to the near-term fate of Eta Car itself and 
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other LBVs — LBVs might fail to lose additional signif- 
icant amou nts of mass and, t hus, never become a Wolf- 
Rayet star (jSmith et al.ll2007l) . 

Only a dozen or so LBVs are k nown in the 
Gala x y ([Humphreys fc Davidson! 119941 : Ivan Gendered 
l2001t IClark et al.1 120051 ). and determining their funda- 
mental parameters and evolutionary status is a key to 
understanding unstable massive stars and Population III 
objects in the high-redshift Universe. 

AG Carinae (=HD 94910, a 2 ooo=10h56mll.6s, 6 2 ooo=- 
60d27minl2.8s) is one of the brighte st and most fa- 
mous LBV stars. iCannon fe Pickering! ((19161 ) classified 
AG Car as a "P Cygni-type" star almost a century 
ago, while i ts pho tometric variability was discovered by 
iGreensteTnl (119381) . I ts high luminosity was first reported 
by iThackeravl ((19501 ) , who also noticed the presence of a 
bipolar nebula around the central star. 

After these pioneering works, intense observing cam- 
paigns in the past decades have provided a wealth of 
observational data on AG Car, showing that the star has 
photometric, spectroscopic, and polari metric variabil- 
ity o n timescales from days to decades dLeitherer et all 
19941: Humphreys fc Davidson] fl99l iStahl et all 120011 : 
van Gendered 120011: iDavies et all 120051 :" see references 
in iGroh et all l2009bf ) . The strongest variability occurs 
cyclically on timescales of decades and is known as the 
S-Dor cycle, during which t he visual magnitude of the 
star decreases to ~ 6.0 mag (jvan Gendered[200l . 

The advent of non-LTE, full line-b lanketed radia- 
tive t ransfer codes, such as CMFGEN ((Hillier fc Milled 
1998), has allowed us to obtain precise fundamental pa- 
rameters for AG Car during the visual minimum epochs 
of the S- Dor cycle and to ga in insights into the nature of 
the star (|Groh et alJl2009bL hereafter Paper I). 

Following the determination of strikingly diff erent ste l- 
lar and wind parameters at different minima ((Paper J ), 
we would like to address the following questions in the 
present paper. Is the bistability mechanism present in 
AG Car? What was the initial mass, and what is the 
current mass of AG Car? Is AG Car close to the Ed- 
dington limit? If so, is it during visual minimum or on 
the way to visual maximum? How does rotation modify 
such a scenario? What is the present evolutionary status 
of AG Car? 

This paper is organized as follows. In Sect. [2] we 
briefly describe CMFGEN, the radiative transfer code 
employed in the analysis of the observed spectra of 
AG Car. A detailed comparison between the minimum 
phases of the S-Dor cycle of AG Car under the light of the 
bistability mechanism are considered in Sect. |3J In Sect. 
|4]we present evidence that AG Car was rotating fast dur- 
ing the visual minimum of 1985-1990. The evolutionary 
status, initial mass, and current mass of AG Car are con- 
sidered in Sect. [5l An analysis of how the fundamental 
parameters of AG Car change from visual minimum to- 
ward maximum of the S-Dor cycle is provided in Sect. [51 
The conclusions of this paper are summarized in Sect. [7l 

2. QUANTITATIVE SPECTROSCOPIC MODELING 

The photometric and spectroscopic observations of AG 
C ar during visual minimum were extensively described 
in iPaperB and references therein, to which we refer the 
reader for further details. 

The detailed spectroscopic analysis of AG Car dur- 



ing visual minimum was done using the radiative trans- 
fer code CM FGEN (|Hilliedll987l [19901 IHillier fc Milled 
119981 11999ft and was presented in IPaper t CMF- 
GEN has been succe ssfully used to stu d y other LBVs 
and relate d objects dSmith et all Il994t iNaiarro et all 
1991 Il997t lNaiajTdl2001t IHillier et al.lll998l I2001L 120061: 
Figer et all 119981: iDrissen et all 120011: iBresolin et all 
2002tlMarcolino et al.ll2007T) . Since the code has been ex- 
te nsively d iscussed in the aforementioned references and 
in iPaperll . we concisely describe its main characteristics 
below. 

CMFGEN assumes a spherical-symmetric, steady- 
state outflow and computes continuum and line forma- 
tion in the non-LTE regime. Each model is defined by the 
hydrostatic stellar radius R*, bolometric luminosity L*, 
mass-loss rate M, wind terminal velocity i>oo , stellar mass 
M, and abundances Z$ of the included species. The ra- 
dius of the hydrostatic core is defined as = i> S onic/3, 
in order to avoid any effects due to the strong wind in 
the determination of R+. The distance to AG Car can- 
not be constrained from the spectroscopic analysis us- 
ing CMFGEN. Throughout this work, we assume that 
AG Car is located at d = 6 kpc dHumphrevs et al]|1989t 
iHoekzema et all 119921 ). We refer to IPaper 11 for further 
discussion on the effects of the distance on the derived 
parameters. 

CMFGEN does not solve the hydrodynamic equations 
to determine the wind structure, and therefore, a velocity 
law v(r) needs to be assumed a priori. In CMFGEN, 
the velocity law is parameterized by a beta-type velocity 
law which is modified at depth to smoothly merge with 
a quasi- hydrostatic structure at the sonic point. The 
velocity structure below the sonic point is iterated to 
fulfill the wind momentum equation with a maximum 
error of 10% for v < 8 kms -1 . Systematic errors in the 
density structure below the sonic point might be present 
due to the proximity of the star to the Eddington limit 
and the choice of M, rotation rate, and viewing angle. 

CMFGEN allows for the presence of clumping within 
the wind. This is accomplished by assuming a homo- 
geneous wind at distances close to the star, and that 
clumps start to be formed at a given velocity v c . The 
wind achieves a volume-filling factor / at large distances, 
as follows: 

f(r) =/+(!-/) exp[-v(r)/v c ] . (1) 

Full line blanketing is included consistently in CMF- 
GEN through the concept of superlevels, which groups 
similar energy levels into one single superlevel to be ac- 
counted for in the statistical equilibrium equations. The 
final atomic model included H, He, C, N, O, Na, Mg, Al, 
Si, Ti, Cr, Mn, Fe, Co, and Ni. Synthetic spectra were 
computed in the observer's frame using CM F_FLUX 
(IHillier fc MiUedH99a iBusche fc Hilliedl2005l) . including 
the effects due to rotation. 

3. THE ROLE OF THE BISTABILITY MECHANISM IN 
AG CAR 

After comparing the lightcurve and the spe ctroscopi c 
evolution during both visual minimum phases ((Paper J ) , 
it can be noted that the star presents different char- 
acteristics in t hese epo chs. An inspection of the AG 
Car lightcurve (IPa per t ) shows that the visual minimum 
phases of 1985-1990 and 2000-2003 are significantly dif- 
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Fig. 1. — Comparison between normalized observed line profiles 
of AG Car obtained during 1985-1990 and 2001 April, around the 
spectral lines of He II 4686 (top panel), He I 5876 (middle), and 
Ha (bottom). The 2001 April observation in the upper panel was 
convolved with a Gaussian profile in order to match the resolution 
of the 1989 March data. The offset in the continuum level from 1.0 
is due to the electron-scattering wings of N II lines at 4630-4653, 
He II 4686, and He I 4713. The data has R ~ 50, 000 in the middle 
and bottom panels and R ~ 3000 in the upper panel. 

ferent, especially regarding their duration. While the 
star was in a visual minimum phase for about 5 years 
between 1985-1990, the duration of the following visual 
minimu m of 2000 -2003 was much shorter, between 2 and 
3 years ([Paper J ). 

There are also significant qualitative differences in the 
spectrum obtained at different visual minimum phases 
of AG Car, which can be inferred by comparing spec- 
tral lines of H I, He I, and He II. Figure [1] displays, for 
each visual minimum, the spectrum obtained during the 
epoch when the maximum value of T e g was achieved, 
corresponding to 1985-1990 and 2000-2001. Figure EQ 
presents the spectra obtained in 1987 June, 1989 March, 
and 2001 April, centered in three key lines for which 
high-resolution spectra were available: Ha, He I 5876, 
and He II 4686. Ha emission, in both the line core and 
the electron-scattering wings, was much stronger in 2001 
than in 1985-1990, indicating that the wind was denser 
in 2001. The P-Cygni absorption component of Ha was 
much deeper in 2001, while it was almost absent in 1985- 
1990. This indicates that the n = 2 energy level popula- 
tion of H was higher in 2001, which was probably caused 
by the lower T e g and, thus, increased amount of neutral 
H. We will come back to this point in Sect. 13.51 The 
He II 4686 emission, on the other hand, was stronger in 
1985-1990 than in 2001, confirming that T e s was higher 
in 1985-1990 than in 2001. The line emission of He I 
5876 was stronger in 2001 than in 1985-1990, but this 
was due to the increased wind density and not because 
of a higher T e g, since the latter would also imply in- 
creased He II 4686 emission, and the opposite is seen. 
The P-Cygni absorption profiles of Ha and He I 5876 
both indicate that the wind terminal velocity was lower 
in 2001 than in 1985-1990. 

The detailed spectroscopic analysis using CMFGEN 
models ([Paper J ) confirms the above scenario and allows 
us to make quantitative conclusions through the compar- 
ison of the physical parameters obtained in the epochs 



when the maximum stellar temperature was achieved in 
consecutive visual minimum phases. While L+ was sim- 
ilar and R* was only 10% different, three other physical 
parame ters were significantly different; namely, T e g, M, 
and Woo ()Paper J ) . Notably, T e g was approximately 5800 
K higher in 1985-1990 than in 2001, while ^ was 2-3 
times higher and M roughly two times smaller in 1985- 
1990 than in 2000-2003. 

Based on our modeling, we suggest that the differences 
in the lightcurve and in the spectrum of AG Car during 
consecutive visual minimum phases are caused by the 
strikingly different wind paramet ers due t o a change in 
the underlying stellar parameters ([Paper J ) . This is con- 
trary to the current ideas about the nature of the S Dor 
cycles, where the visual minimum phase is representative 
of a quiescent phase of the LBV star. In this scenario, 
visual minimum phases should be equal to each other, 
which is not the case for AG Car. Thus, the questions 
to be answered are: Why are the physical parameters 
different in consecutive visual minimum phases? Is there 
any intrinsic reason for this to occur? 

The differences in T c ff, M, and Woo during the min- 
imum phases of AG Car are immediately reminiscent 
of the bistability mechanism seen in line-driven winds 
(IPauldrach fc Puld 11990c ILamers et ail [19951: IvTnk et al.l 
1999), which has already been suggested to be respon- 
sible for the variable m ass-loss rate seen in AG Car 
([Vink fc de Koterl[200l . Before going further and an- 
alyzing whether our results agree with previous works, 
we shall briefly summarize the physics behind the bista- 
bility mechanism. 

3.1. Background on the bistability mechanism 

There are two potential effects related to sudden 
changes in the stellar parameters, which are commonly 
referred to as bistability. 

First, according to the original proposition of 
iPauldrach fc Puli ([1990D when analyzing models for the 
LBV P Cygni, the classical bistability mechanism is char- 
acterized by a change in the optical depth of the Ly- 
man continuum (rL yc ), which in turn changes the ion- 
ization stages of metals in the wind, causing a change in 
the amount of line driving, and thus, ultimately chang- 
ing M and v^. La ter, this concept wa s extended to a 
bistability jump by ILamers et al.l ([1995D to explain the 
physical parameters of blue supergiants, suggesting that 
the winds of those stars have two different regimes of M 
and , which sw itch abruptly f rom o ne to the other 
around 21,000 K. ILamers et al.l (|1995l ) proposed that 
blue supergiants with T c s < 21,000 K could typically 
have M 2 times higher (correspondingly, Voa is 2 times 
lower) than blue supergiants with T e g > 21, 000 K. The 
physical explanation for this phenomenon is a change in 
the ionization str ucture of iron in the inner stellar wind 
(|Vink et al.lll999l ) , which changes the number of spectral 
lines able to absorb/scatter radiation and, thus, able to 
dri ve the wind . This in turn affects the value of M and 
JVink et al l[l999l). 



While the IPauldrach fc Pulsl (I1990D work referred to 
the wind o f a sing le LBV star, ILamers et al.l ([1995D and 
iVink et ail ([l99_9) were analyzing di fferent OB super 



giants . Both ideas were combined by IVink fc de Koterl 
(2002), who suggested that the variable mass loss com- 



4 



Groh, Hillier & Damineli 



monly inferred for LBVs is due to subtle changes in the 
line driving due to recombination/ionization of Fe 3+ into 
Fe 2+ and Fe 2+ into Fe + . However, based on a detailed 
spectroscopic analysis with CMFGEN of a large sample 
of OB stars, the idea of a "jump" in M and Voo of OB su- 
pergiants around T P ff > 21,000 K has been disputed by 
ICrowther et~aT1 ([2006), who suggest a more continuous 
variation of M and Voo with temperature. 

Second, there is also the wind instability, which causes 
sudden changes in the Ha line profiles, intensity of 
Fe II lines, and TL yc , for a relative small change in T*, 
M, and/or T c g. This effect is well known in LBVs 
and has been discu ssed in the lite rature for P Cygn i 
UNaiarro et al.l I1997T). HDE 316285 (IHillier et al.l[l998h . 
and Eta Carinae ((Hillier et al.l 1200 If) . Interestingly, the 
wind instability and the classical bistability appear to oc- 
cur at similar T R ff, causing feedback effects as discussed 
bv lPauldrach fc Pulsl lfl99TO . 

3.2. The behavior of the Lyman continuum of AG Car 

during visual minimum epochs 

In order to check whether the bistability mechanism 
can explain the presence of different stellar and wind pa- 
rameters in consecutive visual minimum phases of AG 
Car, we first turn our attention to the underlying cause 
of the bistability mechanism, which is the change in TL yc . 
We analyzed the behavior of tl yc in selected models ob- 
tained for both minima, which are shown in Fig. O 
While we obtained TL yc < 1 throughout the wind dur- 
ing 1985-1990 (i.e. optically thin), an increase by three 
orders of magnitude in t Ly C was derived in models dur- 
ing 2000-2003, making the Lyman continuum completely 
optically thick. Therefore, we first conclude that the dif- 
ferent stellar parameters at the base of the wind of AG 
Car in different visual minimum phases cause a huge 
change in TL yc , which in turn cause significant differ- 
ences in M and Voo, implying a very different spectral 
morphology. In particular, we conclude that when AG 
Car had T cff > 19,500 K (i.e., during 1985-1990), the 
Lyman continuum was optically thin, M was ~ 2 times 
lower than in 2001, and Uqo ~ 3 times higher than in 
2001, when the Lyman continuum was optically thick 
and T cff ~ 17, 000 K. 

Our results confirm the findings from lVink fe de Koterl 
( 2002) that AG Car lies on diffe rent sides of the bistabi l- 
ity regime at different epochs (jPauldrach fe Pulsl il990). 
We found that between 1985-1990, AG Car had T eff ~ 
22,800 K, which we interpret as being on the hot side 
of the bistability regime, while in 2001 AG Car had 
T ff ~ 17,000 K, which we interpret as the star cross- 
ing the bistability limit into the cool side. This result 
implies that the effective temperature for occurrence of 
the bistability in AG Car is somewhere between 17,000 
K and 22,800 K. Based on our detailed CMFGEN spec- 
troscopic analysis, we suggest that this limit is closer to 
the earlier, around 17,000K - 18,000 K, because in this 
case, small changes of 5% in T*, M, and/or iZ* caused 
significant changes in TL yc , in the ionization structure of 
the wind, and in the emerging spectrum. 

3.3. The iron ionization structure in the inner wind of 

AG Car 
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Fig. 2. — Top: Optical depth of the Lyman continuum (black 
dashed line), Rosseland optical depth (red dotted), and ratio be- 
tween the wind temperature and the grey temperature (blue dot 
dashed) as a function of velocity, derived from CMFGEN models 
of AG Car for the 1985-1990 and 2000-2003 minima. The mod- 
els are ordered by decreasing T B g (and increasing M) from left to 
right. Bottom: Corresponding ionization structure of Fe and H of 
the inner wind of AG Car as a function of velocity. Since H is 
mainly ionized in the inner wind of AG Car at all epochs shown in 
this Figure, only the H° abundance is shown for clarity. For each 
model, the velocity of the sonic point is indicated by a cross at the 
horizontal axis. 



We shall now discuss the iron ionization structure pre- 
dicted by our CMFGEN modeling in the inner wind of 
AG Car, since the huge change in tl vc is expected to 
impa ct the iron ionization structure (|Vink fe de Koter] 
[2001) . 

The inclusion of several physical processes in the cur- 
rent generation of radiative transfer stellar wind codes 
such as CMFGEN turned out to be crucial for the anal- 
ysis of the ionization structure of the wind of early-type 
stars. In particular, the exact treatment of radiative 
transfer in the co-moving frame, the inclusion of full 
line blanketing, wind clumping, and consistent model- 
ing of the wind and photosphere have all had a huge 
effect on the ionization structure of early-type stars (e.g 
Hillier & Mill enfl9p ICrowther et al.l [2002t IHillier et al 



20031: iBouret et al.ll2005t iMartins et al.ll2005H Puls et al 



2006t lHamann et al.l 120061: ICrowther! I2007T T We antici 



pate that the ionization structure of LBV stars will also 
be affected. 

The Fe ionization structure of AG Car during 1985- 
1990 and 2000-2003 obtained from the CMFGEN models 
is presented in the bottom row of Figure [2] We noticed 
that the ionization structure of Fe before the sonic point, 
which is around 15km s _1 for all the above models and 
below which M is set, is modified when t L yc changes. 
During 1985-1990 the dominant ionization stage was 
Fe 3+ , while in the following visual minimum phase of 
2001 there is a steady increase in Fe 2+ as T e ff decreased. 
Since Fe 2+ has a higher opacity than Fe 3+ , particularly 
in the region where the stellar flux is high, more mo- 
mentum can be removed from the radiation field and 
transferred to the wind, causing a higher mass-loss rate. 
Nevertheless, as seen in Fig. [21 Fe 3+ is still the dominant 
ionization stage below the sonic point in the temperature 
regime found for AG Car during visual minimum. 
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3.4. Is there a jump in the wind properties of LBVs due 
to the bistability mechanism? 

iVink fc de Koterl (|2002f ) proposed that the bistability 
mechanism is the cause of the variable behavior of the 
mass-loss of LBVs as a function of the effective tempera- 
ture, predicting the presence of jumps in the M vs. T c ft 
relationship at T e & ~ 21,000 K and T off - 10,000 K 
due to the recombination of Fe 3+ into Fe 2+ and Fe 2+ 
into Fe + , respectively. 

Figure [3] presents M as a function of T e g during the 
visual minimum epochs of AG Car analyzed in Paper I. 
Is there clear evidence of a jump in M as a function of 
T e s as AG Car crosses from one side of the bistability 
regime to the other? From the point of view of the spec- 
troscopic modeling, we cannot give a definitive answer 
to this question for the following reasons. In general, we 
can see a progressive increase in M when T e s decreases 
in the range 21,000 — 14,000 K. Close to the region 
where the Lyman continuum becomes optically thick 
(aka the bistability limit) at T cS ~ 17, 000 - 18, 000 K, 
we can identify AG Car in the two regimes predicted by 
IVink fe de Koterl (|2002f ). On the hot side, M is lower 
than on the cool side, while Woo is higher. However, 
whether the star switches abruptly from one regime to 
another or changes gradually cannot be determined since 
AG Car has never been observed at temperatures very 
close to the bistability limit. Long-term spectroscopic 
monitoring has shown that T c s of AG Car decreases 
quickly from visu al minimum to maximum (on timescales 
of a few months, iLeitherer et~aTlll994l: iStahl et al.ll2001t 
Paper III), and monitoring AG Car in intervals of weeks 
would be required. Even then, our ability to derive the 
stellar and wind parameters would likely be hampered 
by strong time-dependent effects, since the flow timescale 
during these epochs is also of the order of a few months. 
Nevertheless, we can conclude from the short timescale 
that the stellar wind is highly unstable when close to 
the bistability limit, and the inner layers might be ex- 
posed to some sort of fe edback from the unstable wind 
(jPauldrach fc PulslfTOflOh . 

The differences in our results from the predictions of 
IVink fc de Koterl (|2002f) are likely a combination three 
factors: a) the inclusion of full line blanketing and 
wind clumping in our analysis, b) the higher L de- 
rived in our wor k comp ared to the value assumed by 
IVink fc de Koterl (|2002ft . and c) due to the fact that 
the ratio Woo/w CS c is definitely not co nstant when T e g 
changes, contrary to the assumption of IVink fc de Koterl 
(|2002t l. Furthermore, M and v^/v^c of AG Car 
are greatly modified by rapid rotation (|Groh et al.l 
2006), a phenomenon which was not included in the 



Vink fc de Koterl ((20021 ) models. The influence of ro- 
tation will also be variable as a function of Teff, since 
the ratio i w/^rnt, changes fr om visual minimum toward 
maximum (Gro h et al 1120061 : see Sect. SJ). 

The behav ior of v^ , as a function of T e g is more ir- 
regular (see iPaper J ). First, there is a decrease in Voq 
from T cff ~ 23, 000 K until T cS ~ 17, 000 K. Then, a 
steep increase until T e g ~ 16,400 K is seen, and a shal- 
low decrease again until T c g ~ 14,000 K. While the 
overall behavior of might be explained in terms of 
the bistability mechanism through changes in the line 
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Fig. 3. — Mass-loss rate (upper row) and wind terminal veloc- 
ity (bottom row) of AG Car as a function of T c g during visual 
minimum phases. 

driving (|Vink fc de Koterl 120021) . there are a number of 
physical mechanisms which likely have a large impact on 
the behavior of on the cool side: the fast rotation 
derived in AG Car, the presence of a latitude-dependent 
wind and departure from spherical symmetry, the accu- 
racy of the velocity structure near the sonic point, and 
time-dependent effects. 

3.5. The H and Fe ionization structures in the outer 
wind of AG Car and the variable P-Cygni 
absorption profile 

The presence of strong Fe II lines in the spectrum of 
AG Car during the visual minimum of 2000-2003 (which 
were absent during 1985-1990) provide solid evidence 
that the ionization structure of the outer wind of AG Car 
was also different in those minima. Figure [4] presents the 
H and Fe ionization structures in the whole wind of AG 
Car during visual minimum phases as deriv ed from the 
CMFGEN models with the parameters from IPaper 1 It 
can be seen that in addition to the recombination of Fe 3+ 
into Fe 2+ in the inner wind during the 2000-2003 visual 
minimum, there is also recombination of H + into H° and 
of Fe 2+ into Fe + in the outer wind of AG Car. This 
explains the appearance of P-Cygni absorptions during 
the latter visual minimum, which is only possible when 
a significant (10 -3 ) fraction of hydrogen is recombined. 
The same mechanism explains the absence of P-Cygni 
absorption during 1985-1990, since hydrogen was ionized 
thr oughout the whole wind. In addition , as fir st noted 
by IHillier et al.l (fl998l ) and iHillier et al.l (pool , the Fe 
ionization structure in the outer wind is dominated by 
the following charge-exchange reaction with H, 



Fe 2+ + H o- Fc+ + H+ 



(2) 



This reaction implies that hydrogen recombination is also 
responsible for the recombination of Fe 2+ into Fe + , thus 
explaining why there are no strong Fe II lines in the 
1985-1990 spectrum while they are relatively strong in 
data from 2000-2003. 

The marked difference in the ionization structure of the 
outer wind during different visual minimum phases of AG 
Car led us to investigate whether the recombination of 
hydrogen and Fe 2+ into Fe + in the outer wind and the 
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Fig. 4. — Ionization structure of H (upper row) and Fe (bottom 
row) as a function of distance in the wind of AG Car during visual 
minimum ( T c ff decreasing from left to the right) . The region shown 
in Fig. [2] corresponds roughly to < log(r/i? + ) < 0.1. 

consequent appearance of a P-Cyg absorption component 
are also due to the bistability mechanism (i.e. change in 
the optical depth of the Lyman continuum). 

We found that when t Ly C increases, the relative popu- 
lation between the hydrogen levels 2 and 1 also increases, 
and so does the optical depth of the B aimer lines, such 
as Ha and H/3, explaining the increase in the absorption 
of those lines. 

Interestingly, variable P-Cygni absor ption has been ob- 
served in the LB V candidate He 3-519 ([Smith et al.lll994t 
iCrowtherl 119971 ). which might indicate that He 3-519 is 
also close to the bistability limit. Further work is highly 
desirable in order to prove whether this speculation holds 
for LBVs in general. 

4. ROTATION 

The analysis of broad absorptions due to Si IV 4088- 
4116 A in the high-resolution AG Car spectra obtained 
during the 2000-2003 visual minimum has shown di- 
rect evidence, for the first time, that some LBVs ro- 
tate at a high fraction of their critical rotational ve- 
locity ([Groh et al.ll2006D . Using a similar methodology, 
the bona-fide LBV HR Carinae was also sh own to pos- 
sess a rotational velocity close to critical (|Groh et al.l 
I2009aft . Since the Si IV lines were also pr esent in AG 
Car during t he 1985-1990 visual minimum (IStahll 119861 : 
iSmith et al.lll994l iWalborn fe Fitzpatrickll2000D . we an- 
alyzed whether a similar fast rotation could be noticed 
during th ese epochs. Exactl y the same technique as de- 
scribed in I Groh et al.l (|2006[ ) was employed. 

Although the spectra of AG Car obtained in 1989 
March have only moderate resolution around the Si IV 
lines (R~100 kms -1 ), significant evidence of fast rota- 
tion can be inferred from the comparison between the 
observed spectrum and the 2-D models presented in Fig. 
[3 The Si IV line profiles from models without rota- 
tion, convolved to the same spectral resolution as the 
observations, are narrower and have stronger central ab- 
sorption than the observed Si IV line profiles of AG 
Car from 1989 March. This is clear evidence that fast 
rotation is present. The best fit was achieved using 
Urotsini = 220 ± 50 kms" 1 , mainly using Si IV 4088 



as a diagnostic (Fig. [3^,) since Si IV 4116 is severely 
blended with He I 4120 at the data resolution (Fig. [5p). 
We estimate that v rot likely had similar values during 
the whole 1985-1990 visual minimum, since there was 
little va riability in the stellar parameters during these 
epochs ([Paper J ) . Indeed, a tracing of a medium resolu- 
tion p hotographic -plate spectrum of AG Car obtained in 
1985 (|Stahll 119861 ) shows broad Si IV 4088-4116 absorp- 
tions which are very similar to those predicted by our 
rotating models. 

Figure [5b presents the rotational velocity of AG Car as 
a function of i?* for each epoch. Using a least-squares lin- 
ear fit, we obtained v Iot oc i?~ 134±014 , where the quoted 
error only includes the statistical errors. Additional sys- 
tematic errors, which are currently hard to estimate, 
might be present due to model assumptions. A flatter 
dependence of u ro t on i?* might be obtained when includ- 
ing the deformation of the stellar surface from spherical 
symmetry, variation of the stellar and wind parameters 
as a function of latitude, and a more complex treatment 
of the density structure at and below the sonic point in 
the model. As a consequence, we cannot exclude that 
v ro t oc R~ . In addition, our result strongly suggests 
that magnetic fields are not important on the surface of 
AG Car; otherwise, significant deviations from the above 
relationship would be expected (A. Maeder, private com- 
munication). If the magnetic field on the hydrostatic 
surface of AG Car played an important role, one would 
expect that the magnetic field would maintain the surface 
material corotating out to some distance. In this case, 
the surface rotational velocity would behave as that of a 
solid body, i. e., increasing with i?*. This is contradictory 
with our findings. 

The rotational period of AG Car during 1985-1990 
is P rot = 13 ± 2 days, assuming R± = 62 Rq and 
v Iot = 220 kms" 1 . LBV stars have low-amplitude pho- 
tometric variability (the so-called "micro- variability" 
in timescales of d ays to weeks ([Lamers et al.l 1199 
Ivan Gendered [200 ID . Generally, the currently accepted 
mechanism for suc h behavior is non-radial pulsation 
([Lamers et a l.l ll998| ). In the case of AG Car, the micro- 
varia bility period during 1985- 1990 was between 11-14 
days ([van Genderen et all 1 19881 ) . 

The similarity between the values of P ro t and the 
micro-variability period is striking. We suggest that 
they might be related to each other and that the micro- 
variability is modulated by the rotational period of AG 
Car, which would be the case if spots are present on the 
stellar photosphere, for instance. Since the stellar photo- 
sphere of AG Car is extended and located in the stellar 
wind, the presence of spots would suggest that structures 
propagate from the stellar surface into the wind, possibly 
causing effects on line profiles. If the photometric micro- 
variability is indeed modulated by rotation, we anticipate 
that the micro- variability period should be proportional 
to i?~ 2 34 in the case of AG Car, since v ro t is proportional 
to R- 134 for AG Car. 

If rotation is somehow related to the micro- variability 
of LBVs, in principle, their rotational velocities could be 
derived from the micro- variability period once the stellar 
radius i?* is known. This analysis is beyond the scope 
of this paper and is the subject of an ongoing work. It 
is interesting to note that for other LBVs such as HR 
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Fig. 5. — a) Medium-resolution Si IV 4088 absorption profile of 
AG Car observed in 1989 March (solid black line) compared to a set 
of CMFGEN models computed without rotation (solid grey line) 
and with a projected rotational velocity of 150 kms -1 (dot-dashed 
line), 220 kms -1 (dotted line), and 270 kms -1 (dashed line). De- 
spite the medium resolution of ~ 100 kms -1 of the observations, 
we can clearly infer the presence of a high rotational velocity in 
AG Car in this epoch, b) Same as panel a, but for Si IV 4116. c) 
Rotational velocity of the surface of AG Car as a function of i?* in 
log— log scale. 

Car and S Dor, the variations in the micro-variability 
pe riod between visual minimum and maximum is about 
5 (jLamers et al.l [19981) . which is consistent with a rota- 
tional modulation. 

5. EVOLUTIONARY STATUS OF AG CAR 

5.1. Initial mass and location of AG Car in the HR 
diagram 

In this Section we com pare the fundamental parame- 
ters obtain ed for AG Car dPaoer Jl w ith the evolutionary 
tracks from lMevnet fc MaederT(f2003[ ) (hereafter MM03). 
There are a couple of caveats which need to be consid- 
ered: 1) the problematic definition of T e g in LBVs and 
WRs because of the presence of an optically-thick wind 
and 2) the strong, long-term variability seen in AG Car 
during the S-Dor cycle. 

The first issue is a direct consequence of the high mass- 
loss rate of LBVs and WRs, which causes the presence 
of an extended photosphere. The hydrostatic base of 
the wind (i?* in our definition) has a higher tempera- 
ture (T*) than that of the photosphere (Teff, defined at 
tross = 2/3). The relative decrease of T e g compared to 
T* depends on how optically thick the wind is, which is 
a fun ction of M, M, and wind velocity law (e.g.. lLangerl 
1989). Additional complications come from the presence 
of wind clumping and time-dependent effects in LBVs, 
which are b oth not considered by evolutionary models. 

Since the IMM03I evolutionary models do not solve the 
radiative transfer in the extended stellar atmosphere as 
CMFGEN does, they cannot self-consistently compute 
T c g during the LBV and WR phases. Only a rough 



estimate of the effect of the optically-thick wind on T c g 
can be calculated based on t he stellar parameters and on 
the predicted mass-loss rate (lLan gcr 1989; Schalle r~et al.l 
Il992t iMevnet fc Maedenl2005b[ ). Obviously, that result 
will only be approximate and might not be accurate, as 
any changes in the assumed parameters (such as M) will 
bias the comparison between T e g determined from the 
CMFGEN models and T c g estimated by the evolutionary 
models. Therefore, comparing the value of T e ff predicted 
by the evolutionary models with values obtained from the 
atmospheric radiative transfer codes might be misleading 
in the case of stars with dense winds such as LB Vs and 
WRs. In this case, stars with a similar T e s can have very 
different underlying stellar parameters at the base of the 
wind if they have different M, M, and/or Voq. 

The second issue is related to the strong variability 
of AG Car due to the S Dor-type variability. Which 
epoch of the S-Dor cycle should be compared with the 
evolutionary tracks, since the models do not account for 
the S-Dor type variability? The answer to this ques- 
tion is not obvious. Although the current ideas about 
the S-Dor cycle consider the visual minimum phase the 
quies cent phase, followed by eruptio ns during maximum 
(e.g. [Humphr eys fc Davidson! fl994f ) . there is much ev- 
idence which might not corroborate this point. First, 
we detected that VLT in AG Car is higher during visual 
minimum than during maximum (Sect. ^j. Therefore, 
if AG Car is to settle in quiescence during visual min- 
imum, it has to do so extremely close to the Edding- 
ton limit. Second, even if we assume that the star is 
at quiescence during visual minimum, AG Car can have 
different stellar parameters in different minimum phases 
( raper Tj). Which one should be compared with the evo- 
lutionary tracks? Third, some well-documen ted LBVs, 
such as S Doradus (|van Genderen et al.lll997| ) and R127 
(jWalborn et al.1 12008). have been observed at maximum 
for extended time intervals. Fourth, there is histori- 
cal evidence that some LBVs, such as AG Car, have 
presented the S-Dor type instability for over a century 
(|van Genderenll2001[ ). and the key point is how long the 
S-Dor cycles are present. What if the S Dor-type insta- 
bility is present during the majority of the LBV phase 
for these stars? In this case we would just be witnessing 
these unstable, massive stars evolving in real time, and 
they would never really be in quiescence during the LBV 
phase. A definitive answer to this point can only be pro- 
vided when the physical mechanism behind the change of 
the hydrostatic radius of LBVs, which ultimately causes 
the S Dor-type variability, is recognized and included in 
the stellar evolutionary models. 

Keeping these caveats in mind, we present in Figure [S] 
the location of AG Car in the HR diagram during the vi- 
sual minimum phases of 198 5-1990 an d 2000-2001 based 
on the parameters derived in lPaper"H compared to those 
predicted by the evo lutionar y tracks for non-rotating and 
rotating stars from MMQl The He abundance of AG 
Car (Y=0.62. iPaper it ) is a strong constraint on the evo- 
lutionary stage of AG Car, restricting the comparison to 
a very specific part of the evolutionary tracks. According 
to them, the bolometric luminosity of AG Car is consis- 
tent with an initial mass of ~ 100 M@, independent of 
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Fig. 6. — Position of AG Ca r in the HR diagram based on the 
physical parameters derived in Paper] for the visual minimum of 
1985-1990 (x s ymbol) a nd 2000-2003 (open triangle). Evolution- 
ary tracks from MM03 for non-rotating stars (a and c) and stars 
with an initial rotational velocity of 300 kms -1 (6 and d) are shown 
for initial masses of 60 Mq (full line), 85 Mq (dotted), and 120 Mq 
(dashed) . The upper panels ( a and b) present the value of TV from 
the CMFG EN models with the effective temperature predicted by 
the IMM03I not corrected by optical-depth effects, while the lower 
panels ( c and d) compar e the value of T e g from CMFGEN with the 
values from the IMM031 evolutionary models corrected for optical- 
depth effects due to the presence of a dense wind. The thick part 
of each track corresponds to epochs where the helium abundance 
is 0.5 < Y < 0.7, which is the range commonly found in LBVs. 



the initial rotational velocit 



Technically, only the|MM03 non-rotating models in the 
range 40 < M < 60 M (25 < M < 45 M Q for rotat- 
ing models) have a bona-fide LBV phase; in models with 
M > 60 Mq the star ski ps the LBV phase and goes 
directly to the WR phase (|MM03f ) . While both the ro- 
tating and non-rotating models can reproduce fairly well 
the effective temperature of AG Car (Fig. [6l lower pan- 
els), only the non-rotating model is consistent with the 
value of T* derived for AG Car (upper panels). The 
evolutionary models with rotation predict a much hot- 
ter ( T* ~ 49, 000 K) WN star at the base of the wind , 
corresponding to a WN7 spectral type (|Crowtherll2007D . 
Even though the stellar evolution model with rotation 
can fairly reproduce the T e g of AG Car, it cannot re- 
produce the value of that we inferred based on the 
spectroscopic analysis. Therefore, based on current evo- 
lutionary models, we can rule out that the progenitor of 
AG Car had a high initial rotational velocity. 

The non-rotating models simultaneously reproduce the 
values of T+, T e g, and He/H; the value of T* is consistent 
with a WN11 spectral type. However, the fast rotation 
of AG Car is puzzling in the non-rotating scenario. One 
possibility is that the initial rotational velocity of AG 
Car was significantly lower than 300 kms^ 1 . Another 
possibility is that additional processes need to be taken 
into account by the evolutionary models in order to fully 
explain the stellar properties during the LBV phase. For 
instance, it would be desirable to include giant outbursts 
of several Mq in the evolutionary models to investigate 

4 Rot ating models are slightly more luminous than non-rotating 
models (MM03T) , but not sufficiently enough in order to change this 
conclusion. 



5.2.1. 



5.2. Current mass 
The minimum mass according to the Eddmgton limit 

The influence of the stellar mass M and effective grav- 
ity g on the spectrum of OB supergiants is usually ob- 
tained from the broadening of H and He lines. However, 
the photosphere is extended when a dense stellar wind 
is present and, as the location of the photosphere ap- 
proaches the sonic point, line emission from the wind fills 
out the original photospheric absorption profile. Thus, in 
the case of LBVs and WR stars, the dense wind causes 
the photosphere to be located at velocities larger than 
the sonic point and strong emission lines appear in the 
spectrum. Therefore, the influence of log g through the 
analysis of hydrogen absorption lines is not applicable to 
AG Car. 

To analyze the influence of log g on the spectrum of AG 
Car, we ran CMFGEN models assuming different stellar 
masses in the range 40 Mq < M < 150 Mq, keeping 
the other model paramete rs constant, such as a He mass 
fraction of 0.62 (jPaper J ). The stellar mass affects the 
Eddington parameter T since 



T = kL*/(4itGcM) 



(3) 



where n is the total wind opacity and the other sym- 
bols have their usual meaning. For epochs where the 
wind density was large enough that w p hot > v s and the 
star was not so close to the Eddington limit (r < 0.7), 
such as during 2002 March - 2002 July and 2003 Jan- 
uary, we found little influence of log g (and hence M) on 
the emerging spectrum. The physical parameters derived 
from models with different log g were virtually the same, 
and little insight on the value of M could be obtained. 
This result is similar t o what has been found for Eta Car 
(Hilli er fc Millerlll998l) and WR stars in general. 

However, for epochs when T > 0.8 (1985-1990 
and 2000-2001), we found that the emerging spectrum 
changed significantly as a function of the adopted M, 
and so did the parameters obtained from the spectro- 
scopic analysis. Unfortunately, M could not be directly 
derived from the models, since a reasonable fit of the 
spectrum could still be obtained by adjusting other pa- 
rameters such as T e ff, M, and /. We derived that the 
changes in M affected mostly T c ff , which varied between 
19600-23800 K when M varied from 40 to 150 Mq, as 
shown in Fig. [7] The value of T e g is sensitive to the 
density and velocity structure below the sonic point, and 
we determined that T c ff is usually 1000-2000 K higher 
in models which adopted a fixed isothermal scale height 
to match a hydrostatic struc ture (using a procedure de- 
scribed by Hillier et al.| [2003). The dependence of T e ff on 
M in this parameter regime can be understood as being 
due to the formation of a progressively more extended 
photosphere as the star approaches the Eddington limit. 
When M decreased (while L+ and T c g were kept con- 
stant) we found that higher ionization was achieved in 
the wind. Therefore, a decrease in T c g is needed to fit 
the observed spectrum again. 

The middle panel of Fig. [7] shows T, obtained from 
our CMFGEN modeling, as a function of M, computed in 
two locations of the wind: at high optical depth (t = 100 
and v ~ 0.03kms _1 ) and at R+ (v ~ 5kms _1 ). Ac- 



On the nature of AG Car II. Evolution close to the Eddington and bistability limits 



9 



cording to our CMFGEN models, AG Car has to have 
M > 40 Mq to ensure T < 1. Different relationships 



were obtained: while T oc M 



-0.96±0.09 



at v ~ 5km s 



a much less steep variation was obtained at the base of 
the wind, at v ~ 0.03km s^ 1 (T oc M-°- 37±0 - 02 ). The 
different steepness at different depths in the wind reflect 
the different photospheric extensions and opacities as a 
function of M. 

5.2.2. Effects of fast rotation on the minimum mass of 
AG Car: the QT limit 

The high rotationa l velocity derived for AG Car dur- 
ing visual minimum ([Groh et al.l 120061 : Sect. [3} will re- 
duce th e maximum luminosity tha t can be reached before 
T = I (|Maeder fc Mevnetl l2000al hereafter MM00) and 
will require a larger value of M in order for the star to 
be below the Eddington parameter modified by rotation 
(To). Here we show that, making a few reasonable as- 
sumptions, we are able to obtain a lower limit for Tq and 
M when takin g rotatio n into account. 

According to IMMOOI . 



k(0)Z*[1 + C(0)] 



AttcGM 



(l 0—^ 



(4) 



where 6 is the latitude angle and ( is a correction factor to 
the von Zeipel gravity-darke ning law to take into account 
the baroclinicity of the star ( MMOOf ) . 

Since the stellar parameters of AG Car were obtained 
in iPaper II through a one-dimensional radiative transfer 
analysis, and the high rotational velocity will cause devi- 
ations from spherical symmetry, there are a few caveats 
which need to be taken into account to obtain a reason- 
able value of Tq . 

First, the mass-absorption coefficient and radiative flux 
may vary as a function of latitude in a fast-rotating 
star. Unfortunately, given the 1-D nature of our radia- 
tive transfer analysis, we cannot infer directly how these 
changes will impact the estimation of Tq. To obtain a 
rough estimate of how close AG Car is to the Eddington 
limit modified by rotation (the so-called fir-limit), we 
ignore latitudinal variations of the mass- absorption co- 
efficient and radiative flux in Eq. |4l which then reduces 
to 

* : T n . (5) 



f 



2irGp„ 



To properly include realistic latitudinal variations of the 
opacity and radiative flux, one would need detailed hy- 
drodynamical modeling coupled with radiative transfer 
calculations in 2-D, which is certainly well beyond the 
scope of this paper. 

At first glance, the assumption that the derived cor- 
responds to the equatorial value of might imply that 
we are neglecting gravity-darkening, i.e, that the flux of 
a rotatin g star is larger at the polar direction than at the 
equator ()von ZeipeJ Il924f ). However, the value of of 
a rotating star derived using a 1-D model also depends 
on th e viewing ang le dCollins 1963, 1965; IGillich et al.l 
2008). For instance, for a Main-Sequence O-star with 
T e g = 35,000 K rotating at 0.9 of its critical velocity 
for break-up, the actual value of is about 25% higher 
than the derived value of L* when observing equator-on 



and using a 1-D model (Warren & Hillier, private com- 
munication) . The presence of the wind will likely change 
this dependence. Takin g into account th at AG Car is 
se en nearl y equator-on (jGroh et al.l 120061 ). obtained 
in IPaperll is likely lower limit, and it is reasonable to as- 
sume that the flux at low latitudes is not much different 
from the one predicted by a spherica l symmetric model 
with L± — 1.5 x 10 6 L Q as in IPaper 1 . 

Equat ion [5] can be rewritten as a function of u rot as 
(IMMOOft : 

r o — r E^r-r > ( 6 ) 



1 



R 2 



where R c is the equatorial radius at a given w rot , i? e b and 
i?pb are respectively the equatorial and polar radius at 
the classical break-up, u c rit,i = (GM/Reh) ' 5 is the first 
critical velocity for rotation, and V is the volum e of the 
star relative to that of a sphere with i? p b, as in MM00. 

R 2 ra 

Assuming that the factor V-0- is close to 1 in AG CaiQ 

and using the definition of w C rit,i, Eq. breads as 



(7) 



1 



2 ".ot-Rpb 

3 GM 



The second caveat regards how R+, determined from 
the 1-D radiative transfer model, relates to i? p b- We 
argue that i?* ~ R p b because 1) w r ot/w C rit,i is at maxi- 
mum 0.7 for AG Car, meaning that the shape of the star 
does not significantly departs from spherical symmetry, 
and 2) the effects of fast rotation on the stellar struc- 
ture were not included in CMFGEN in the computation 
of the hydrostatic atmospheric structure below the sonic 
point. Using this simplification, we can then replace i? p b 
by R* and, using T and v mt determined from the spec- 
troscopic analysis and assuming i — 90°, obtain Tn from 
models with different M using Eq. [7J The results from 
this calculations are shown in the bottom panel of Figure 

m 

Since v vot is considerable, the minimum mass of AG 
Car will be larger than the value of ~ 40 Mq derived in 
Section 15 . 2 . II according to the T-limit. Based on the 1-D 
CMFGEN models we conclude that, ultimately, the fast 
rotation of AG Car requires M > 60 Mq in order to sat- 
isfy the Eddington limit modified by rotation. It would 
be desirable to confirm this result with appropriate full 
2-D, line-blanketed, non-LTE radiative transfer models 
when these become available. 

5.3. The evolutionary mass budget and current mass of 
AG Car: a massive star close to the Eddington 
limit 

The first constraint to the current mass of AG Car 
comes from evolutionary model s. Acco rding to the non- 
rotating evolutionary models of lMM03l . and assuming an 
initial mass of 1 00 Mq an d the derived surface He abun- 
dance (Y=0.62, IPaper J ), the current mass of AG Car 
would be roughly 60 Mq. This value is extremely close 



5 MM00 show that this factor varies between 1 and 0.813, and 
our test calculations show this produces a negligible effect com- 
pared to the other errors associated with the analysis. 
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Fig. 7. — Top: Derived value of T e g as a function of the stel- 
lar mass for different CMFGEN models for the 1985-1990 visual 
minimum phase. Middle: Eddington parameter T as a function of 
mass for different CMFGEN models for the 1985-1990 visual min- 
imum phase. Note that the T e ff and M of each model were also 
changed in order to fit the spectrum. The different symbols corre- 
spond to values of T obtained at high optical depth (r = 100 and 
v ~ 0.03km s _1 (asterisks) and at v ~ 5km s — 1 (open circles). A 
least-squares linear fit is shown for each sample. Bottom: Edding- 
ton parameter modified by rotation Fq as a function of the stellar 
M computed at the base of the wind at v ~ 0.03km s — 1 (open 
squares). The arrows pointing upwards indicate that QV > 1, 
which was obtained for M < 60Mq . 

to the minimum mass of AG Car according to the mod- 
ified Eddington limit derived above, which would imply 
that AG Car is at the flT limit. 

There are two important caveats of doi ng such a com- 
parison. First, the evolutionary models of MM031 do not 
include giant eruptions, while AG Car is known to have 
had a strong mass ejectio n ~ 10 years ago. Second, 
the evolutionary models of MMOli include the mass- loss 
recipe during the Main S equence (MS) an d Blue Super- 
giant (BSC) p hases from (IVink et al.ll2001l ). while recent 
findings (e.g., iBouret et alj|2003t see Sect. [T]) seem to 
suggest a lower mass-loss rate during these epochs. If 
that is true, the evolutionary models would predict that 
AG Car would reach the LBV phase with a lower M. 

A second constraint to the current M can be obtained 
by computing the mass budget throughout the evolution 
of AG Car. If we take into account that AG Car has 
a massive ne bula (15 — 30 Jlfc) composed of material 
from the star ([Lamers et al.l[l998l ) and an initial mass of 
100 M Q , then the mass budget is tight: AG Car would 
have M ^ 70 — 85 M© assuming zero mass loss during 
the MS and BSG phases. A more realistic estimate of 
M - 2 - 4 x 10~ 6 M Q yr~ 1 during — 3 x 10 6 years would 
lead to M ~ 60 — 65 M© . If the mass-loss rate during the 
MS and BSG phases was even lower, we could make the 
case for AG Car having M ~ 70 M Q . Therefore, evo- 
lutionary considerations and the presence of a massive 
nebula suggest that AG Car has M ~ 60 — 70 Mq and is, 
thus, dangerously close to the Eddington limit modified 
by rotation. 

Obviously, a precise determination of the amount of 
nebular mass around AG Car and detailed evolutionary 
models that include a lower mass-loss rate during the MS 
than the one currently used, and including the mass lost 



in brief giant eruptions, are highly desirable. 

6. ON THE MODUS OPERANDI OF THE S-DOR CYCLE 
TOWARD MAXIMUM 

The evolution of LBVs during their S-Dor cycles 
is commonly assumed to be due to changes in the 
stella r radius at consta n t lum i nosity dLeitherer et al.l 
1989t iHumphrevs et all 119891: iLeitherer et al.l 119941: 



Smith et al. 1994; Shore et al. 1996; de Koter et al. 



1996 



Stahl et all 120011: IVink fc de Koterl 12002). which imply 
changes in the effective temperature. Based on a de- 
tailed spectroscopic analysis of multi-epoch observations 
of AG Car using CMFGEN, we argue that this scenario 
needs some revisions in order to explain the tem poral be - 
havior of the physical parameters of AG Car (jPaper J ) 
which are shown in Fig. |SJ 

We found that the stellar radius of AG Car changes 
during the S-Dor cycles, increasing from visual mini- 
mum to maxi mum, as has been determined by previ- 
ous works (e.g iStahl et alJuOOlD . However, unlike pre- 
vious works, using CMFGEN we can precisely derive 
both the hydrostatic radius of AG Car defined as 
v — fsonic/3) and the photospheric radius (i? p hot, defined 
at tr oss = 2/3) for a given epoch, and thus, compare how 
they change as a function of time. We found that both 
R± and i? p hot increase from visual minimum to maxi- 
mum, ruling out that the S-Dor type variability is caused 
only by an expanding pseudo-photosphere due to an in- 
crease in M (Lcithcrcr et al.l 119851 : lLamers fc de Loorel 
[19871 IDavidsonl fl987l) . The increase in M is not high 
enough to explain the changes in i? p hot! definitely, our 
results show that a physical expansion of i?* is needed 
to explain the S-Dor type variability from visual mini- 
mum to maximum in the lightcurve, as first suggested 
by Ivan Gendered (|1982f) . The relatively high M derived 
for AG Car compared to B supergiants in the same tem- 
perature range, in combination with the proximity of the 
star to the Eddington limit, causes the unavoidable for- 
mation of an extended pseudo-photosphere, which trans- 
lates into different val ues for R * and iZphot- 

A major result from lPaper J is that the bolometric lu- 
minosity of AG Car decreases during the expansion of 
the hydrostatic radius, which is contrary to the afore- 
mentioned status quo explanation for the S-Dor variabil- 
ity. The dete c tion o f a decrease in is not unique to AG 
Car; lLamers! (|1995f ) also detected that another famous, 
highly-variable LBV, S Doradus, presen ts a sim il ar be- 
havior. We support the hypothesis by lLamers! (|1995l ) 
that the decrease in L± is due to the energy needed to 
expand the outer layers of the star. Vari ations in L+ 
were also detected in the LBV AFGL 2298 (|Clark et al.l 
l2009h . However, in this latter case, the variation of L+ 
were apparently not related to changes in i?*. 

Another striking result is that the Eddington param- 
eter r of AG Car decreases as the star expands and its 
and T e ff decrease (Fig. [5]). Such a reduction is, at 
least in part, caused by the decrease in L+. For AG Car, 
r ~ 0.78 during the hottest epochs, and decreases to 
r ~ 0.70 as the star evolves toward maximum. 

The Eddington parameter modified by rotation flT is 
even more dramatically reduced from visual minimum to 
maximum (Fig. [S]) , since rotation will be relatively more 
i mportant when the star has a higher T e g (see discussion 
in lGroh et al. 2006). While AG Car is significantly close 
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Fig. 8. — Stellar and wind parameters of AG Car ijPaper 111 as 
a function of the stellar radius obtained in different epochs as the 
star moves toward maximum. From top to bottom, we show the 
evolution of (a) stellar temperature at the base of the wind (T*) 
and effective temperature (T e g), (b) bolometric luminosity (L*), 
(c) photospheric radius (-Rphot)i (d) Eddington parameter (T), (e) 
Eddington parameter modified by rotation (f2r), (f) projected ro- 
tational velocity (i> ro t), (g) mass-loss rate in units of f0 — 5 Mqjt -1 
(M), and (h) wind terminal velocity (■Uoo)- 

to the Eddington limit modified by rotation (fir ~ 0.93 
during 1985-1990 and 2001 April) during visual mini- 
mum, we found that flT ~ 0.72 as the star expands. 
Therefore, if the proximity to the Eddington limit re- 
flects the instability of the outer layers of AG Car, the 
star is more unstable during visual minimum than during 
maximum. 

The evolution of the wind parameters of AG Car is 
regulated by the bistability mechanism (Sect. [3]). How- 
ever, the behavior of the m ass-loss r ates and wind ter- 
minal velocities obtai ned in iPaper II do not c losely fol- 
low the predictions of lVink fe de Koterl (|2002ft . This is 
due to several reasons: our modeling determined that 
AG Car is twice as m assive as the value derived by 
IVink fc de Koterl (I2002T) . more luminous, and closer to 
the Eddington limit. Furthermore, AG Car has a high 
rotational v elocity which vari es from visual minimum to 
maximum (IGroh et all [20061). Since the models from 
IVink fc de Koterl (|2002t l do not include rotation, signifi- 
cant differences are anticipated. 

7. CONCLUSIONS 

The detailed spectroscopic analysis of AG Car using 
the radiative transfer code CMFGEN has provided a 
wealth of information on the LBV phenomenon and on 
the S-Dor type variability. Below, we summarize the 
main conclusions of this paper. 

1. The high rot ational velocity of AG Car detected by 
IGroh et all (|2006) during the visual minimum of 2001- 
2003 was also obtained during the previous visual mini- 
mum of 1985-1990, amounting to v rot sini ~ 220km s _1 . 



We found that the rotational velocity is proportional 
to ^- 1 - 34 ±°- 09 ) but including systematical errors which 
might arise due to the star being extremely close to the 
Eddington limit, we cannot rule out that v Iot oc R+. 

2. We noticed that the rotational period of AG Car 
during the visual minimum of 1985-1990 (P ro t = 13 ± 2 
days) is remarkably similar to the cyclical photometric 
micro- variability of 11- 14 days observed dur i ng th ese vi- 
sual minimum epochs ()van Genderen et~aLl fl988h . We 
suggest that the micro-variability period of AG Car 
might be modulated by the stellar rotation. If rotation 
is also responsible for the micro- variability seen in other 
LBVs, in principle, their rotational velocities can be de- 
rived using this me thod. 

3. We found in IPaper II that the consecutive visual 
minimum phases of 1985-1990 and 2000-2003 are dif- 
ferent in duration, maximum magnitude achieved, max- 
imum stellar temperature, mass-loss rate, and wind ter- 
minal velocity. We suggest that these differences arise 
due to different stellar parameters, which cause different 
mass-loss rates and wind terminal velocities. We explain 
this behavior in terms of the bistability mechanism of 
line-driven winds, which is ultimately regulated by the 
optical depth of the Lyman continuum. From 1985-1990, 
AG Car was on the hot side of the bistability ( T^g =22800 
K), while the star was on the cool side (T e ff=17000 K) 
from 2000-2003. Therefore, AG Car is the first star ever 
confirmed to be on different sides of the bistability at dif- 
ferent epochs based on detailed spectroscopic analysis. 

4. We noticed that AG Car has never been observed 
exactly at the boundary of the bistability limit, suggest- 
ing that significant instability is present when that hap- 
pens. Although we identify AG Car in two regimes, the 
absence of observations close to the boundary of crossing 
from one side to the other of the bistability does not allow 
us to confirm the presence of a jump in the wind prop- 
erties of AG Car. An increase in M as a function of T e g 
was found as the star crosses the bistability limit, while 
Uoo decreases and presents a more complex behavior. 

5. The duration of the last two visual minimum phases 
of AG Car and the variability inside the visual mini- 
mum are related to the maximum temperature achieved 
at each minimum and, therefore, on which side of the 
bistability the star was. It seems that visual minimum 
phases that have T c ff high enough to put the star on the 
hot side of the bistability limit have a longer duration, 
while visual minimum phases which have T c g low enough 
to put the star on the cool side, are brief. This behavior 
needs to be verified by further spectroscopic and photo- 
metric monitoring of AG Car during the next decades. 

6. We suggest that the proximity to the Eddington 
limit modified by rotation (flT) plays a key role in ex- 
plaining the behavior of AG Car. The high value ob- 
tained for its luminosity (L± = 1.5 x 1O 6 L0), together 
with the hi g h rota tional velocity of AG Car detected by 
IGroh et"al~l (j2006[ ) and by this work, put the star ex- 
tremely cl ose to the Eddington limi t modified by rota- 
tion, QT (|Maeder fc Mevnetl l2000al ) . For an assumed 
mass of 70 M , flT ~ 0.93 during 1985-1990 and 2001, 
and decreases to ~ 0.73 as the star expands towards max- 
imum. 

7. Despite the extreme conditions found in AG Car due 
to its high luminosity and fast rotation, the mass-loss 
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rate remains relatively modest and does not wildly in- 
crease due to the proximity to the Eddington limit. Sur- 
prisingly, the mass-loss rate is even quite modest when 
compared to Eta Carinae, which has M two orders of 
magnitudes higher than AG Car, but a luminosity that 
is only 3 times higher. 

8. Due to the proximity of AG Car to the Eddington 
limit and the mass loss, an extended pseudo-photosphere 
is formed, hiding the hydrostatic surface of the star. 
Even if the mass-loss rate increases towards maximum, 
that is not sufficient to explain the changes in the spec- 
trum. An expansion in the hydrostatic radius is required 
in order to reduce the effective temperature. 

9. Contrary to the current paradigm, the maximum 
value of the Eddington parameter is found during vi- 
sual minimum, not when the star is evolving towards 
maximum. Therefore, the picture of a quiescent star at 
visual minimum followed by an eruption during visual 



maximum due to an increase in the Eddington parame- 
ter might not correctly describe the S-Dor cycles. 
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